Electric polarization in correlated insulators 
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We derive a formula for the electric polarization of interacting insulators, expressed in terms of the full 
Green's and vertex functions of the system. We use the formula to investigate changes in the electric polariza- 
tion of the half-filled ionic Hubbard model in the DMFT framework. Interaction renormalizes both ionicity, and 
covalency of the system and competition between these renormalized quantities determines the electric polar- 
ization of a correlated band insulator. Based on a quasi-particle approximation, we comment that the LDA+U 
based calculations, which does not account for the quasi-particle weight renormalization, will overestimate the 
electronic part of the electric polarization. This is sizable when the quasi-particle weight is far from one for 
active bands in the polarization. 
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INTRODUCTION 

The electric polarization, P, is a measure of the position 
differences between the center of the band electrons and the 
lattice ions. It is generally non zero for a material which lacks 
inversion symmetry and plays the role of the order parameter 
in the theory of ferroelectricity. Over the past two decades, 
there have been important conceptual and computational ad- 
vances in the first principles calculation of this quantity. [1] 

The modern theory of polarization expresses P in terms of 
the Berry phase acquired by a Slater determinant, describing 
the insulating state as an effective one particle theory. [2] A 
complementary picture of the polarization, as the displace- 
ment of the Wannier centers of the single particle states in- 
volve, emerges naturally in this formalism. [3] 

The Berry phase formalism was extended beyond an effec- 
tive single particle picture, by considering the changes of the 
phase of the many body function of all the electrons in the 
solid as a response to the changes in the boundary conditions. 
[4, 5] The resulting expression can also be recast in terms of 
the resolvent of the Hamiltonian of all the interacting electrons 
in the many body ground state. [6] The many body wave func- 
tion formalism to calculate polarization was applied to simple 
model Hamiltonians, [7] but is prohibitively difficult to carry 
out in practice for realistic models of the electronic structure 
of a solid. 

In this paper we formulate the problem of the calculation 
of the electronic polarization of a correlated electronic system 
in terms of the one particle Greens function and vertex func- 
tions. We derive a general expression for the polarization, and 
discuss its implementation within dynamical mean field the- 
ory (DMFT). [8] We apply the formalism to a simple model of 
a correlated insulator. The approach has similarities to recent 
studies of incorporating correlations in the calculation of the 
topological indices of topological insulators [9, 10], and to the 
early work of Volovik [11] and gives useful insights into the 
effects of Hubbard correlations on the electric polarization. 

In an extended system the change in polarization, AP, in- 
duced by distorting an atomic coordinate, Ri Q — > Rj Q + 
£ARi Q , that breaks inversion symmetry is given by the in- 
tegrated bulk transient current as the system adiabatically 



evolves from the initial state (£ = 0) to the final state (£ = 1), 
i.e., [1] 



AP= / d^= [ d£-L ( drJ(r,0 
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(1) 



where J(r, £) denotes the current density and V^ e ii is the 
primitive-cell volume. Implicit in the analysis is that the sys- 
tem must remain insulating everywhere along the path, as oth- 
erwise the adiabatic condition fails. 

Assume the system is described by H (£) for a given value 
of the distortion. A small change in £, i.e., £ — ^ £ + <5£, per- 
turbs the Hamiltonian of the system by (SH/6£)S£ = H'^SS, 
and induces a current. In order to evaluate the current explic- 
itly, we assume S£(t) = <5£ exp(— iut+rjt) with rj an infinites- 
imal positive number and finally we consider the v — > limit. 

According to the Kubo formula the current density induced 

Clir^fa = Q, V )S& where 



by a perturbation is 3(v) 

Cjff' (q = 0, z/) is q = component of the Fourier trans- 
form (in space and time) of the retarded correlation function, 
[12] 



a 



j(r) ff '(r') = -mt-t')([k(r,t),H' ( (r', f)]>f. (2) 

Here the subscript £ for operators shows that they are in the 
Heisenberg representation with respect to the instantaneous 
Hamiltonian and (■•■)$ means that the average is taken with 
respect to the instantaneous spectrum. In the limit of v = 
(static distortion), cjf^q = 0,v = 0) = 0, therefore the 
driving force is <5£ and J oc 6£. We evaluate J as 

J(q = 0,0 =i(|-q j ( | ) «(q = 0,^)1^ 

where we used analytic continuation in the last line and we de- 
fined transport coefficient C^, . Substituting Eq. (3) in Eq. (1) 
gives an expression for AP. 
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Note that ^tC^Jj, , is anti-symmetric under exchange of j and 

H' operators, i.e., = — 5ft£^j, which can be seen from 

its Lehmann representation (see supplemental material). [13] 
In a diagrammatic series expansion we can express Cy^f 1 
in terms of interacting Green's functions and irreducible 
particle-hole vertex functions. We can find the perturbation 
expansion for C^h 1 by differentiating with respect to iv n . This 
introduces a derivative vertex {d/diuj m ) (see supplemental 
material). [13] The diagrams contributing at the C$h> are of 
two types: (i) diagrams which are separated into three pieces, 
each of which has a vertex, by cutting three electron lines. We 
denote the contribution of these diagrams to the polarization 
by APi. (ii) All the rest diagrams which their contribution to 
the polarization is denoted by AP2. The change in polariza- 
tion from the diagrams of type (i) is given by 

ap 1= - r dt(—)YY 

Tr ( A.f G«)A$G«)^^G«) - \f o a|4) 

where we have used the identity {dG^/diu> m ) = 
—G^(dG^~ 1 /diu} m )G^ s> and the anti-symmetric prop- 
erty of CjH> in writing Eq. (4) down. Aj, Ah> are the particle- 
hole vertex functions and f3 is the inverse temperature. The 
interacting single particle Green's function of the system is 

G ( £(iu m ) = [(iw„+/z)l-H^(k)-£^(k> m )]-S (5) 

where Hq (k) denotes the non-interacting part of the Hamil- 
tonian. While and are gauge covariant, AP is 

gauge invariant. However, the select of the interaction matrix 
and the expression of the current do depend on the gauge. [ 14] 
Here we write Hq in the periodic gauge. [15] si (k, iui m ) 
denotes the electron self-energy. Bold quantities are n x n 
matrices where n denotes number of orbitals within the unit 
cell. 

One can use the Ward identity to replace the electric cur- 
rent vertex by Aj (e) (k) = -(dG^liiWmj^/dk) and the dis- 
tortion current vertex by A^|)(k) = — (9G|^ {iujm)^ 1 / d^). 
[16] Therefore, the Eq. (4) can be rewritten as 

ap 1= - f 1 (—) y y -e^- 

*(^G<0^G<0^G«)) (6) 

\ Ok Ot; OlUJ m J 

where we have introduced an anti-symmetric tensor, £ k ? w ™. 
Eq. (6) follows from Eq. (4) due to the anti-symmetric prop- 
erty of the derivative relative to k and £ and cyclic property 
of the trace. In the non-interacting case Eq. (6) reduces to the 
Berry phase formula (see supplemental material). [2, 13] 



It is not possible to write AP2 in the form of Eq. (4). How- 
ever, we have evaluated it perturbatively in various models and 
found it to be very small. 

In the following we will work in the DMFT approxima- 
tion, where we ignore the current vertex corrections while 
keeping the self-energy corrections associated with elec- 
tron correlation. In DMFT, AP2 vanishes and the corre- 
lation function reduces to the bubble contribution with the 
fully interacting Green's function and the following ver- 
tices; (dG^-i/dk) = -(dH^/dk), (dG^^/d^ = 
-(dU^/dO - (dS^/dO, and (dG^-^diLOn) = 1 - 
(0Efc>/0iw n ). 

To investigate the influence of Hubbard correlations on the 
electronic part of the polarization and to benchmark our for- 
malism we turn to a simple model of ferroelectric materials, 
the ionic Hubbard model (IHM) in one dimension. On a bi- 
partite lattice the IHM Hamiltonian is defined by the following 
equation [17-20]: 

H = - 1 5}1 + (-m](clc l+la + H.c.) 

io 

+ A^(-l)' i n i(T + C/^n lt n 4 (7) 

ia i 

where, t(> 0) and U denote the hopping amplitude between 
nearest-neighbor sites and the on-site Coulomb repulsion, re- 
spectively. — A(+ A) denotes the local potential energy for the 
A(B) sublattice, respectively. 

The model is relevant for the study of organic ferroelectric s 
with stacks of alternating donor and acceptor molecules. The 
paraelectric phase is centrosymmetric (so that all A — B dis- 
tances are the same), whereas in the ferroelectric phase there 
is a relative displacement of the A and B sublattices, which 
breaks the inversion symmetry and produces an alternating 
pattern of short and long A — B distances. Here we modify 
the hopping amplitude by ±£ alternatively to include a dimer 
term that breaks the inversion symmetry. 

At half filling, the noninteracting undistorted system (£ = 
0) is a band insulator (BI) with a charge gap equal to 5 C = 2 A. 
Upon turning on and increasing U the charge gap, 6 C , in the 
correlated band insulator shrinks until a discontinuous transi- 
tion to a Mott insulator (MI) occurs with a hysteresis region, 
if the anti-ferromagnetic (AF) long-range order is not allowed 
to set in. In the Mott phase a large gap is established. [21, 22] 
As the system is driven deeper into the Mott phase, the gap 
size increases. Both the BI and the MI are non-polar at £ = 0; 
however, depending on how one defines the zero reference 
point of the polarization one can assign a polarization value to 
these states. Here we define the zero reference point to be the 
usual band insulator where both electrons occupy the Wan- 
nier function centered on the lower energy site, then using the 
classical point charge model we assign P^/ea = -1/2 to a 
Mott insulator with no long range order. In the centrosymmet- 
ric system without spontaneous inversion symmetry breaking 
[24], increasing U from zero causes a charge exchange be- 
tween sites. However, charge flows symmetrically, the macro- 
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FIG. 1. (Color online) Electric polarization as a function of distor- 
tion, £, for several value of U in the PM phase (For U/t < 7.5 
the system is in the band insulator phase with zero reference point 
of polarization P(£ = 0)/ea = 0. For U/t = 8.0 data is 
shown in the Mott phase with zero reference point of polarization 
P(£ = 0) feet = -1/2). Black dashed lines show Eq. (13). 



scopic current is zero and therefore AP' cl ) = 0. In the pres- 
ence of the AF order, system has a transition to the AF phase 
at a smaller critical U. [23] 

Next we present our results in the presence of a distortion. 
We set A = 0.5t and use Eq. (6) and DMFT to evaluate the 
electronic part of the electric polarization of the distorted IHM 
in the paramagnetic (PM) and AF phases. The Bloch repre- 
sentation of the non-interacting part of the IHM Hamiltonian, 
Eq. (7), is given by 

H (fc) = -2tcos(ka/2)T x + 2£tsm(ka/2)T y - At z (8) 

where T XtVtZ are the Pauli matrices and a is unit cell length. In 
the DMFT approximation the current vertex and the distortion 
vertex are given by 

= at sm(ka/2)T x — a£tcos(ka/2)T y , (9) 



dk 



-2t s\\\(ka/2)r y — 



(10) 



where X is a diagonal matrix with elements Y,a and Eg. 

Fig. 1 shows the polarization of the system, P^ cl y(£)/ea, 
obtained from an evaluation of Eq. (6) for several values of U 
as a function of £. As can be seen in the correlated BI (U < 
U c ~ 7.5*), the absolute value of the P$j(£)/ ea increases 
with £ and saturates at —1/2. The saturation value happens at 
smaller £ upon increasing U (note that for small values of U 
the saturation value is not apparent on this figure, as it occurs 
at larger £). In the non-interacting case, a system with small 
ionicity is more polarizable. In order to explain results with 
non-zero U with U < U c , one needs to understand how the 
correlation renormalizes the ionicity and the covalency of the 
system. To get some insight we work in the QP approximation 
and derive an analytical equation for the electric polarization 
at small dimerization. The QP Green's function is given by 

G«)(k,a;) = zV>„l - H$(k)]-V/ 2 , (11) 



where H$(k) = z x /2( H ^(k) + 5tX[f(0) - i ul)z 1 /2 i s 
the QP Hamiltonian. The quasi-particle weight is denoted by 
z ~ [1 — 3S(icL>o)/wo] _1 , where luq = ir//3 is first fermionic 
Matsubara frequency. Substituting Eq. (11) into Eq. (6) and 
using the fact that dG~^/diuj n = z~\ one can see that an 
extension of Eq. (SI 3) in the supplemental material [13] to 
the interacting IHM is obtained by replacing bare quantities 
with renormalized ones, A — > A™ and t —> zt. In our case 
z is a diagonal matrix with za = zb = z. The renormal- 
ized ionicity can be obtained from the renormalized charge 
gap which in turn can be obtained by identifying the poles of 
the renormalized propagator, Eq. (11), and is given by 



S qp 



o a ren 



2z[A - 3f?(E A (iu;o) - Ss(iw ))]. (12) 



The frequency range of the validity of the Fermi liquid as- 
sumptions decreases upon increasing U, but the charge gap 
also shrinks simultaneously. A closer investigation of the self- 
energy and local density states of the correlated insulator show 
that one can continue to assume the Fermi liquid assumptions 
on the frequency range between the highest pole of the valance 
band and the lowest pole of the conduction band for all inter- 
action strengths. Therefore, Eq. (11) describes low energy 
physics precisely and the Eq. (12) is the true charge excitation 
gap, S c = Sf. 

Thus, in the QP approximation the polarization for small 
distortions is given by: 



P 



(el) 



A,U 



(Ar7(^))sin 2 (fca/2) 



N 



^ [(A™"/(zt)) 2 +4cos 2 (fca/2)] 3 / 2S 

(13) 

which shows that for small £ the electric polarization of the in- 
teracting system is the same as a non-interacting system with 
effective ionicity given by the charge gap of the system renor- 
malized by quasi-particle weight, (A 1 q '/ ) n /(zt)). Upon increas- 
ing U, (A^/t) decreases faster than z which leads to a re- 
duction of (Aqp n /(zt)) which in turn leads to a more polar- 
izable system where a small distortion can trigger substantial 
changes in the electric polarization. In Fig. 1 the QP result 
Eq. (13) is shown with dashed lines. The agreement between 
the QP approximation and the full answer is very good in the 
linear regime of polarization. Notice that at strong interaction 
strengths Eq. (13) fails because it is built on the assumptions 
that 5RE is linear and 5E vanishes identically on the whole 
bandwidth. 

Fig. 2 shows P&ir(£ = 0.08i)/ea for several values 
of U obtained from the full- and QP-calculations. In the 
QP-calculation we employ H^p to define an effective non- 
interacting system and calculate P using the Berry phase for- 
mula. [2] It is worth mentioning that same results for polariza- 
tion can be obtained by working with (H^ (k) + 3?s[f (0) - 
fil) as the QP Hamiltonian, which provides the correct ratio 
of the renormalized ionicity to the covalency. However, the 
latter Hamiltonian gives incorrect results on the charge gap 
and does not account for quasi-particle weight renormaliza- 
tion. Consistently with Eq. (13), the QP results agree very 
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FIG. 2. (Color online) Electric polarization for a distorted system 
with £ = 0.08t for several interaction strengths, U, using the full 
calculation and the quasi-particle (QP) approximation. 



well with the full calculations at weak to intermediate interac- 
tion strengths. However, at strong coupling, P obtained from 
the QP method deviates from the full answer because it does 
not account for the bandwidth renormalization correctly. Inset 
of the Fig. 2 shows the charge gap as a function of U. 

At large interaction strength, for example U = 8.0t, the 
system is in a Mott phase. In a Mott Phase, the variation 
in the electronic part of the polarization as a function of the 
distortion is very small (here we restrict ourself to small £). 
Indeed in a Mott phase, charge fluctuations are strongly sup- 
pressed and an ionic displacement induces only a very small 
current. In the Mott phase the effective non-interacting sys- 
tem described by H 9P is an extremely ionic system where one 
sublattice is fully occupied and the other one is empty. In re- 
ality, however, in the true ground state the charge is almost 
uniformly distributed in the system. Thus, the QP approxima- 
tion fails in the Mott phase. 

Next we investigate the AF phase which allows us to com- 
pare our results with the QMC results on IHM. Investiga- 
tion of the staggered magnetization, m = (jia-\ — n ai) = 
— (fiBf — ns\), of the centrosymmetric structure, £ = 0, 
as a function of the interaction strength shows that at zero 
temperature the system shows long-range AF for U > 2 At. 
The onset of the AF phase coincides with the onset of the 
possible spontaneous bond order phase. [25] Fig. 3 shows 

\ P A%^ = °- 02 ) - p a%(£ = °)l obtained from an im p!e- 
mentation of Eq. (6) in the AF phase along with a QMC cal- 
culation reproduced from [7]. In the QMC data, the phase dia- 
gram of the undistorted ionic Hubbard model includes a spon- 
taneous dimerized phase with non-zero polarization which 
means that the reference point of polarization, y(£ = 0)> 
has different value than us. However, the shown data is AP 
which is independent of the reference point. Also, the QMC 
calculation was done for a slightly different A. Nevertheless, 
the two calculations agree very well with one another. 

In conclusion, we have introduced a practical many body 
approach to the calculation of the electric polarization of in- 
teracting insulators. We have implemented and tested the for- 
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FIG. 3. (Color online) [Pj^(£ = 0.02) - P%$(£ = 0)j obtained 
from an implementation of Eq. (6) in the AF phase along with a QMC 
calculation reproduced from [7]. Inset shows a schematic phase dia- 
gram of the system. 



malism in the context of a model Hamiltonian and shown that 
while the electronic polarization is affected by Hubbard cor- 
relations, it is a less sensitive quantity well approximated by 
a quasi-particle approximation and depending on the ratio of 
the gap and the covalency of the quasi-particle band struc- 
ture. These results justifying the success of the traditional 
electronic structure of methods which do not include an ex- 
plicit frequency dependence of the self energy . Notice how- 
ever, that in some correlated materials such as HoMnC<3 this 
method has been reported to overestimate the electric polar- 
ization. [26, 27] This can be understood as resulting from the 
effects of the bandwidth renormalization on the polarization 
described here. The formalism of this paper can be imple- 
mented in realistic LDA+DMFT codes. This together with 
a comparison with accurate experimental ARPES studies of 
correlated ferroelectrics are interesting open problems. 

We are grateful to F. Marsiglio, M. E. Pezzoli, A. 
Soluyanov, G. Ortiz, K. Haule and D. Vanderbilt for useful 
discussions. This work has been supported by NSF DMR- 
0906943. 
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Electric polarization in a correlated insulator - Supplemental materials 



I. ANTI-SYMMETRIC PROPERTY 



The anti-symmetric property of the transport coefficient appearing in Eq. (3) can be seen from the Lehmann representation of 
the correlation function. In the general case we have for the correlation function [1] 

Z ^— ' V + £j n — tj m + IT] 
nm 

where Z = Tr[exp(— 0H)]. Therefore, Cab = K^^ab)\ u —q nas following form 

r 4 1 (n\A\m){m\B\n) , a E -/3E m \ 

Now if we commute A and B operators, we get 



_ _ • 1 ( n \ B \ m )( m \A\n) u _ pEn _ o _p E . 

Comparing Eq. (S2) and Eq. (S3), gives 



r aA - ,1 \- {n\U\m)(m\A\n) , » En 0Em) 

CbA ~ l Z^(E n -E m + ir l )^ e 6 > (S3) 



KCab = -KCba (S4) 
SCab = ^Cba (S5) 

since only real part is giving the current, then one conclude that the transport coefficient is anti-symmetric. 



II. DIAGRAMMATIC EXPANSION OF THE CORRELATION FUNCTION AND ITS DERIVATIVE 

Fig. SI shows diagrammatic series expansion of the C-fj, . Lines show the full Green's function and r(kio; m + 
iv n , piuj m r, piuj m i + iv n ,\iiu] m i) is irreducible particle-hole interaction vertex. Bare vertexes are given by Aj (k) = 
(9Hq /9k) and A]|;(k) = (9Hq respectively. By inserting a (d/div n ) vertex and then looking at the v n — > limit we 
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can find the perturbation expansion of the (9Cjjj, /diis n ) | _ Q (see the second and third row of the Fig. SI). Aj and Ah' are 
vertex functions and the last row shows their diagrammatic expansion. 

The diagrams contributing at the (dC-jj, (q = 0, iv n ) / div n )\ Vn —Q can be classified into the following two types: (i) diagrams 
which are separated into three pieces, each of which has a vertex, by cutting three electron lines, (ii) diagrams which cannot be 
separated into three pieces by the same operation (See the third row of the Fig. SI). 

(k, iui m + iu n ) (k, iu m + iv n ) (p, ioj m , + iv n ) (k, iu m + iv n ) (p, iw m i + iu n ) 

C«',(q = 0,^„)= Aj< - 



(k,iw m ) (k,«w„,) (p,«%') (k,iw m ) (p,ji 

<T^ x ' : < TFT> x Aj <TFT!> A " + Aj <TTH"> 



<::> 

o 



(i) (ii) 



FIG. SI. (Color online) Diagrammatic expansion of the correlation function, C^,{(\ = 0,iv n ), and its derivative, (9Cj^,(q = 
0, iv n )/divn) \u n =o- Lines show the full Green's function, T is irreducible particle-hole interaction vertex and T' = (dT / div n )\ v , 
where v n is bosonic frequency. Dash line shows d/diuj n vertex. Aj and A H i are the particle-hole vertex functions and the last row represents 
their Bethe-Salpeter equation. 



III. NON-INTERACTING SYSTEM 



Here we show that in the non-interacting case the Eq. (4) reduces to the Berry phase formula. Using the band representation 
of the Greens function, g^^(ioj m , £) = [iuj m l — e^J] -1 , and 9i Wm g^(iw m , = 1 one can rewrite Eq. (4) as 



- -^EE^(8£(^0(^)8S(M„,e)(^)gS(^e) 

kcr w m 

ie \"\^\"\" V Mm\ d£ I g£ Mini dk j™ 



EEEE 



ku n m w m " e ikrj)(* a; ™ ~ e mkcr)(* w ™ — e ik<r) 

« r 3H S) r aH ^) -( 9H £) r SH *o 



^ kcr n m c,e „kcr V^m e nkirA !w m e mk<W 



By summing over Matsubara frequency and using the fact that the Fermi function is one for occupied bands and zero for 
unoccupied bands, we find (terms with n = m vanish) 



(IP^ C ' ie v-^ \ -« x - { gk" )nm\ jmn ~ ( g^—)nm\ g^" )mn rets 

~dT = ~ N ^ ^ ^ 73® _AO )2 (S7) 

kcr nEoccm^tn V c rikcr c mkcr/ 



which reduces to the Berry curvature [2] 
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= -f E E [^A iVk«2L> - v k (^i^ii>] (ss) 

ktr nGocc 

^ing Al^l«£L> = (eit - eSj<V k «®>«L) and (^tl^l^L) = (4L ~ ^L)(^4tl^L) 
n 7^ to. These expressions can be obtained by taking the k(£)-gradient of H^|ic? ) = e^to-l^nko-) an( ^ taking the inner 
product with J. Here, IVvoL-) = elk r l u lkcr) aTe tne Bloch states and the summation is performed over the occupied 
bands. The right hand side of Eq. (S8) is the Berry curvature which is gauge invariant and thus observable. [2] The Berry 
curvature is nonzero in a wide range of materials, in particular, in crystals with broken time-reversal or inversion symmetry. In 
the periodic gauge [3], the integral on £ can be done analytically and it leads to a two-point formula for the electric polarization 
that involves only the initial and final states of the system [4]: 



A p(el) =p( el )(f / )-P(«*)(6), 
N 



p (el) (o = f E E(4iiv k ki)- <S9) 



kcr n£occ 

Ionic Hubbard model In 1D-IHM, the energy bands are e±(k) = ±e(k) = ±[A 2 + 4i 2 cos 2 (£;a/2) + 4£ 2 i 2 sin 2 (fca/2)] 1 / 2 
and the bare current vertex and the bare distortion vertex in the band representation are given by (dH^(k)/dk) = 

V- 1 (k)[-atsm(ka/2)T x + a£t cos(te/2)r y ]U(fc) and (dfl$\k)/dg) = U _1 (fc)[2t sin(/sa/2)x 2/ ]U(fc), respectively, r are 
Pauli matrices. U(fc) and U~ 1 (k) diagonalize the IHM Hamiltonian as Ho(fc) = U(fc)e(fc)U~ 1 (k) and are given by 

H (fc) = [-2tcos(ka/2)r x + 2^tsm(ka/2)r y - At z ], (S10) 

U(fc) = 1 [2^tsin(fco/2)] + (e(k) - A)r x + 2t cos(fca/2)r z ], (SI 1) 

v /2e(fc)(e(fc) - A) 

U _1 (fc) = 1 [-2ffisin(fca/2)1 + (e(fc) - A)r x + 2<cos(fca/2)r z ]. (S12) 

^2e(k)(e(k) - A) 

At small £ the denominator of Eq. (S7) is given by (e_ k(T — ^\ a ) 2 — 4(A 2 + 4i 2 cos 2 (fca/2)) while the imaginary part of 

the numerator is /dk) 12 {dH^ } /d£) 21 - {dH { ° /^)i 2 (5H^ /dk) 2 i 2t 2 sm 2 {ka/2) • (A/^A 2 + At 2 cos(fco/2)). 

Therefore, the polarization at small £ is given by 

p(«Qftn~ 2ea V (AA)sin 2 (fca/2) 
*A,ffW- at [( A /i)2 + 4cos 2 (fc a /2)] 3/2 
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